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I. INTRODUCTION 



The study of the energy levels of simple atomic systems (muonium, positronium, hydrogen 
^ \ atom, muonic hydrogen and others) with high precision plays significant role for the check of 
' the Standard Model and the bound state theory with very high accuracy. The two-particle 
bound states represent important tool for the exactitude the values of fundamental physical 
constants (the fine structure constant, the electron and muon masses, the proton charge 
radius etc.) The observation of thin effects in low energy physics of simple atoms 

can be considered as necessary supplement to the construction of large particle colliders 
for deep penetration to the structure of elementary particles and search of new fundamental 
interactions. Such atomic experiments can improve our knowledge about elementary particle 
interactions on small distances what may be reached only at very high energies |2j]. 

The effects of strong interactions play essential role in the energy spectrum of the muonic 
hydrogen just as electronic hydrogen. On one hand, they are connected with two electromag- 
netic proton form factors (electric Ge and magnetic Gm) describing the distributions of the 
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electric charge and magnetic moment. In the Lamb shift case the main contribution to the 
energy spectrum of order iZa)^ is determined by the proton charge radius which appeares 
as differential parameter of these distributions. So, the comparison of the experimental data 
and theoretical value for the Lamb shift obtained with the corrections of high order over a 
gives the effective approach to obtain more reliable value of the r^. The measurement of the 
IP — IS Lamb shift in /xp with the precision 30 ppm allows to obtain the value of the proton 
charge radius which is an order of the magnitude better in the comparison with different 
methods The calculation of the nuclear structure corrections in the hyperfine splitting 
of the energy levels (see 0, IE & IE S llll) can be done only on the basis of whole 



proton electromagnetic form factors. Last experimental measurements of the form factors 
Ge and Gm were carried out in Mainz 20 years ago fl^. 

On the other one, important contribution of stron g in teractions to the hydrogen spec- 
trum is connected with the proton polarizability [l3[ ll4L Il5l . Ilfil |. It appeares already in 
the one-loop amplitudes of the muon (electron) proton electromagnetic interaction when 
different baryonic resonances can be produced in the intermediate states as a result of the 
virtual Compton scattering on the proton. Exact calculation of such effect can be done by 
means of experimental data and theoretical models for the polarized nucleon structure func- 
tions. The proton structure and polarizability effects lead to main theoretical uncertainty in 
the expressions for different energy levels including the hyperfine splitting of the hydrogen 
ground state: 



3 (mi+m2)3' 

where fip is the proton magnetic moment in nuclear magnetons, mi is the muon mass, 1712 
is the proton mass, S'^^^ represents the QED contribution, S^^^ is the contribution of 
hadronic vacuum polarization (HVP), the corrections S^^^ and 5^°' are the proton struc- 
ture and polarizability contributions. The expression (1) is valid both for the muonic and 
electronic hydrogen but the exact value of these corrections is essentially different for such 
atoms. The ground state hydrogen hyperfine splitting measurement was made many years 
ago with very high accuracy |l7| : 

Aug^'^iep) = 1 420 405.751 766 7(9) kHz. (2) 

Existing difference between the theory and experiment without accounting the proton po- 
larizability contribution can be expressed as follows [l8| : 

Ai?S(e p) - AE^^^ie p) ,,,,,, , 



E^ie p) 



-4.5(L1) X 10~^ (3) 



This quantity contains one of the main uncertainties connected with inaccuracies of the 
proton form factor determination. Dominant part of the one-loopproton structure correction 
is defined by the following expression (the Zemach correction) : 



AEz = E 



p 2/ia f dp 



[p2 + 1^2)2 



Ge{-p')Gm{ -P 



2^ 



E^{-2i^a)Rp, W = afi, (4) 



where /i is the reduced mass of two particles, Rp is the Zemach radius. In the coordinate 
representation the Zemach correction (4) is determined by the contraction of the charge 
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Pe and magnetic moment pm distributions. The Zemach radius represents the integral 
characteristic of the proton structure effects in the hyperfine sphtting of the energy levels. It 
may be considered as new fundamental proton parameter in the hydrogen atom. Numerical 
value of the Zemach contribution is equal 

/^Ez = -1.362 ± 0.068 meV, (5) 

where the 5% estimation of the uncertainty is connected with the measurement of the proton 
electromagnetic form factors So, the measurement of the muonic hydrogen hyperfine 

splitting as for the electronic hydrogen with similar accuracy 30 ppm as in the case of the 
Lamb shift can give new information about possible value of the contributions 5'^*'" and 5^"' 

Hi. 



Such experiment demands corresponding theoretical study of different order corrections 
with the same precision. Analytical calculation of the hydrogen hyperfine splitting was 
carried out during many years [18, JOy and reached the accuracy 10~*. But these calculations 
can not be used directly for the muonic hydrogen after the replacement the electron mass 
to the muon mass. The reason consists in the proton structure effects. Indeed in the case 
of the muonic hydrogen the dominant region of intermediate loop momenta is of order the 
muon mass. So, the calculation of higher order amphtudes with good accuracy can be 
based only on their direct integration with the account of experimental data on the proton 
electromagnetic form factors. 

The investigation of different contributions to the energy levels of the muonic atoms was 
done many years ago in Ref. 0|. So, at present there is need for new more complete 
analysis of all possible corrections in the HFS of the fip with the declared accuracy 30 ppm. 
Main corrections of order to the hyperfine splitting of the 25* state in the pp were studied 
m Ref.Q. They are very important for the extraction of the Lamb shift value 2P — 2S in 
the experiment. In this study we calculate different contributions of orders and to the 
muonic hydrogen HFS which are determined by the effects of electromagnetic and strong 
interactions. The aim of the work consists in obtaining the numerical value of the ground 
state HFS in the muonic hydrogen with designated accuracy which can serve as reliable 
guide for corresponding experiment. Some basic problems of the HFS measurement in the 
muonic hydrogen were discussed in Ref. [i^ . 



II. EFFECTS OF VACUUM POLARIZATION IN THE ONE-PHOTON 

INTERACTION 

Our calculation of different energy levels of the hydrogen-like atoms are carried out on 
the basis of the quasipotential a ppr oach where the two-particle bound state is described by 



the Schroedinger-type equation |24| : 



where 



El-ml = El-ml 



Pr = E1E2/M is the relativistic reduced mass, M = Ei + E2 is the bound state mass. 
The quasipotential of the equation (6) is constructed in the quantum electrodynamics by 
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FIG. 1: Effects of the one- and two-loop vacuum polarization in the one-photon interaction. 



the perturbative series using projected on positive states the two-particle off mass shell 
scattering amplitude T at zero relative energies of the particles: 

y = V^(i) + V(2) + v(3) + ...^ T = TW+T(2)+r(3) + ..., (7) 

We take the ordinary Coulomb potential as initial approximation for the quasipotential 
Vip, q, M): Vip, q, M) = V^{p- q) + AV{p, g, M). 

The increase of the lepton mass when we change the electronic hydrogen to the muonic 
hydrogen leads to the decrease of the Bohr radius in the fip. As a result the electron Compton 
wave length and the Bohr radius are of the same order: 

n 

— : = 0.737384 

fie'' rrieC 

(me is the electron mass, /z is the reduced mass in the atom fip). An important consequence 
of last relation is the increase the role of the electron vacuum polarization effects in the 
energy spectrum of the fip 2^|. The effects of the vacuum polarization in the one-photon 



interaction are shown in Fig.l. 

To obtain the contribution of the diagram (a) Fig.l (the electron vacuum polarization) 
to the interaction operator there is need to make the following substitution in the photon 
propagator j2^: 



1 a /-i ^'(1-t) 

72 ^ - / dv ,,,^ ^ (9) 
P TT Jo k^{l - v^) - Ami 

At (— /c^) = fc^ ~ iil{ZaY ~ ml^ZaY (electronic hydrogen, /Xg is the reduced mass in 
hydrogen atom) we obtain — a/lSvrmg omitting first term in the denominator of right part 
of Eq.(9). But when k'^ ~ /i^(Za)^ ~ m\{ZaY (muonic hydrogen, mi is the muon mass) 
than /ia and m^ are of the same order and it is impossible to use expansion over a in the 
denominator of Eq.(9). To construct the hyperfine part of the quasipotential in this case 
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(the muonic hydrogen) in the one-photon interaction we must use exact expression (9). We 
take into account that the appearance of the electron mass mg in the denominator of the 
amphtude leads effectively to the decrease the order of the correction. It is well known that 
the hyperfine splitting quasipotential has the form j26l |: 



yHFS^^^ = ^l±^-i[(.,.,)k^ - (.,k)(.,k)]. (10) 

' mim2 4 

For the S'-states 



(11) 



rHFS/u\ _ SttZq: (Tl(T 2, 

3mim2 4 

/«=1. 792847337(29) is the proton anomalous magnetic moment. Averaging the potential 
(11) over the Coulomb wave functions we obtain main contribution of order [Za)'^ to the 
HFS of the IS'-state in the system fip (the Fermi energy): 

= -(Z a)^-^ (1 + k) = 182 AASmeV. (12) 
3 mim2 

The modification of the Coulomb potential due to the vacuum polarization (VP) is de- 
termined by means of Eq.(9) in the momentum representation as follows (25l |: 

.c„.>_ " r (2f + 1) 



In the coordinate representation we obtain: 



VUr) ^ ^ r .; V^5 + ^ f-£^e-- ) . (14) 



37r Ji 

The contribution of the electron vacuum polarization to the hyperfine splitting part of the 
I7 quasipotential for the S'-states can be derived in a similar way in the momentum and 
coordinate representations: 



Vi^, yp(k) = — — -(crio-2)k - / ao.2\ ^^^ (1^) 



^^^"^^ " 3mim2 3^^ ^ 



vr5(r) - 



(16) 



Using Eq.(16) we can obtain the electron vacuum polarization correction of order to 
the HFS in the iip. Taking the wave function of the IS'-state 

^ioo(r) = ^e-^^ W = fiZa, (17) 
we represent this correction in the form: 



AE,,ypHFS= " \ ^ ^ ' --^ / ^ + l (18) 



t2 

-2 
.2^ 



6 



X 



1 



0.374 meV. 



The contribution of the muon vacuum polarization (MVP) can be found by means (16) after 
the substitution nie mi. This correction is of order due to the reason mentioned above. 
Numerical value is equal 



AE^/ljyp = E^—-^Za^ = 0.002 meV. 
'' lb mi 



(19) 



The diagrams of the two- loop electron vacuum polarization shown in Fig.l (b,c,d) give the 
contributions of the same order a^. The interaction operator corresponding to the loop after 
loop amplitude can be obtained using the relation (9). In the coordinate representation 



^l-y, VP-VPV ) 



SnZall + k) {(Ti(T2) ( a 



3mim2 
V^f^{2r]^ + 1) 



-dr] 



5{v) 



m 



T:r{rf — ^^) 



(rj 



X 



(20) 



and the contribution to the energy spectrum 



AE^^^^yp_yp = 0.001 meV. 



(21) 



To calculate the contributions of the diagrams b, c in Fig.l which are determined by the 
polarization operator of the second order it is necessary to make the substitution in the 
photon propagator ^27.] : 



1 



a 



TT 



2 ^1 



fiv) 



Ami + ^^(1 



a\2 2 /-i 



dv= -] - 



dv 



IT J 3 Jo 4mg + fc^(l— f^) 



X 



(22) 



Li'; 



x|(3-w2)(l + ^;2 



1 + v 



+ 2Li2 



1+v 



+ - In — In — — In — 

2 1-v 2 1 



Inf 



+ 



In 



+ 



-v(3 — f ^) In 

2 ^ ' 4 



To find numerical value of this correction we write the quasipotential in the coordinate 
space: 



/\yHFS \ 

^^l7, 2-loop VP) 



iirZail + k) (a 



3m\m2 



TT 



2 fl 



f{v)dv 

(1-^2) 



5(r) 



mz 



7rr(l — f 2) 



The potential (23) gives the contribution to the HPS in the muonic hydrogen 

^Eii,%ioop VP = 0.002 meV. 



(23) 



(24) 



We calculate all contribution numerically and the results are presented with the accuracy 
0.001 meV. 
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III. SECOND ORDER OF THE PERTURBATION THEORY 



The corrections of the second order of the perturbative series in the energy spectrum are 
defined by the reduced Coulomb Green function (RCGF) |28'l]: 

Gi(r,r') =Y.~g^i{Ty)YUr^)Y:^{W). (25) 

l,m 

The radial wave function gni{f'i r') was obtained in Ref. 0| as an expansion over the Laguerre 
polynomials. For the - state 

,,o(r,r) = -4;.Za(^E^^^^^^^-^ + -----je ^, (26) 

where x = 2fiZar, are the Laguerre poljTiomials: 

Some terms of the quasipotential contain the 6{f) so we have to know the quantity Gi{r, 0). 
The expression for the RCGF was found in this case in Ref . jl^ on the basis of the Hoestler 
representation for the Coulomb Green function after the subtraction the pole term: 

Gis{r,0) = -^-1 \2x{\nx + C)+x' -bx-2] , (28) 

where C = 0.5772... is the Euler constant. The main contribution of order in the second 
order of the perturbation theory can be written in general form: 



oo 



^-^1 SOFT - 2^ ' V^^J 

n=2 ^1 ~ 

where /SV^^^ ~ 5(r). Using the relations (14), (28) we can present Eq.(29) as follows: 



x/ (ixe"''^^+^U2x(lnx + C) + - 5x - 2 



0.734 meV. 



The contribution of order in the second order of the perturbative series which is de- 
termined by the vacuum polarization can be derived from Eq.(29) changing AV-^^^ — *• 
AV^^yp. Using exact expressions for the wave function ^/'^(r) (17) and the RCGF (28) we 
write this correction 



<^ l+2Fri^x (31) 



OO 
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a b 

FIG. 2: Proton structure corrections of order [Zaf' . Bold circle in the diagram represents the 
proton vertex operator. 



mi 



me 
W 



rj 



(1 + 



m^X \ ^ / 1 _|_ rriErj ^ ^ 
W \^ ^ W 



W I W_ _|_ 



In 



1 



2(1 + ^ 



(1+ 



'W 



y l-^ 



Vmeg me?? 

1 

1^ 



W 

me»7 



+ 



(32) 



1 + 



meg 



+ 



1 



w 



1 - In 1 + 



W 



+ 



Numerical value of this contribution is equal 



^E^soPT = 0-002 meV. 



(33) 



The second order of the perturbative series gives also other relativistic corrections of order 
(Za)^ including recoil effects which were studied in Ref. |j3fli,i31|,,3^]. Corresponding numerical 
data are in the Table 1. 



IV. PROTON STRUCTURE AND VACUUM POLARIZATION EFFECTS 

The proton structure corrections in the system fip are relatively large in the comparison 
with the electronic hydrogen. In the HFS of the muonic hydrogen these corrections are 
defined in the leading order by the one-loop diagrams in Fig.2. 

To construct the quasipotential corresponding to these diagrams we write the proton 
tensor: 



P2 - k + m2 
{P2 — ky — 1712 + ^0 



7.^1 - 7^(r,xk''F2 
2mo 



u{p2 



(34) 
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where P2i(l2 are four momenta of the proton in initial and final states. The construction 
of the potential can be essentially simplified using the projection operators for the system 
muon-proton on the states with definite spin: 



75, 7r(^^i) = [u{p2)v{pi)]s=i 



:i+7°) 

2^2 



(35) 



where is the polarization vector of the state with the spin 1. Neglecting relative motion 
momenta of the particles in the initial and final states we obtain 
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32k^ 
ml 



+lQk%Fl + l2^k%FiF2 + QAk^klF^ + ?,2k^Fl + 6Ak^FiF2 
Transforming the integration in Eq.(36) to the Euclidean space 

d'^k = An I k^dk / sin^ (j)d(j)^ ko = k cos 



{k^ - Amlk^ik^ - Amlk^oY 



(37) 



we make analytical integration over the angle and present the correction (36) as one 
dimensional integral over the variable k: 



AE 



HFS 
str 



-E' 



Za 



^7m^(l + K 



-6, 



10 



(38) 



V{k) 



2Fie 



fx 



mim2 _ m2)k{k + J Anil + ^ 



USF^ml - l2SFiF2m{ + lQF(k!'+ 



2i2 



2 32F^m^A;^ AF'ik'' AF^k^ 



+6AFiF2k' + 16F^k^ + 



ml 



+ 



2 

mf 



ml 



+ 



■X 



X 



(mi — m2)k{k + yAml + k'^) 
128F^ml + 128FiF2ml - IGF^k^ - 64^1^2^^ - ASF^k^^ 



To cancel infrared divergence in Eq.(38) it is necessary to add the contribution of the iteration 
term of the quasipotential (10) in the HFS of the fip: 



A ^ffps rr ff P <? 6A u'^ (Z ttf (l + k) f^dk 

AEZ%^. = - < V^i^ X X Vr, >fj'= -- ^ ^ ^ 3 / I^' 

3 mim2iin-^ Jo A; 



(39) 



where the angular brackets represent averaging of the interaction operator over the bound 
state Coulomb wave functions and the index HFS shows the hyperfine part in the interaction 
term of the quasipotential (10). The sum of the expressions (38) and (39) coincides with 
the result of Ref. i^]. The integration in Eqs.(38) and (39) was done by means of the 
parameterization of the proton electromagnetic form factors obtained from the analysis of 
elastic lepton-nucleon scattering 2| . Numerically the proton structure correction of order 
(Za)^ is equal 

AEf,f' + AE^l',, = -1.215 meV (40) 
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abed 

FIG. 3: Vacuum polarization and proton structure corrections of order a{Za)^. Dashed line in the 
diagram represents the Coulomb photon. 



Moreover, the effects of the proton structure must be taken into account carefully in the 
amplitudes of higher order over a shown in Fig. 3. 

The contributions of the diagrams (a) and (b) in Fig. 3 to the potential can be found as 
for the amplitudes in Fig.2. taking into account the transformation of one exchange photon 
propagator as in Eq.(9). Corresponding correction to the HFS of the energy level is equal 



A T^HFS 



Za 



1 



l-^]dv 



87r(l + K)n3 TT Jo P(l-t;2)+4m2 



dkVvpik), 



(41) 



where the potential Vvp{k) differs from V{k) in the relation (38) only by the factor k'^. 
Despite of the finiteness of the integral (41) the amplitude terms of the quasipotential in Fig. 3 
(a), (b) must be completed by two iteration terms shown in Fig. 3 (c), (d). First addendum 
< X G-l" X AVyp^ > of order a(Za)^ must be subtracted because the 27 amplitudes (a) 
and (b) in Fig. 3 produce lower order contribution. Second term < Vyp x x V^^^ > 
which is also of order a;(Za)^ has the structure similar to Eq.(29) of the second order of 
the perturbative series. The contributions of discussed iteration terms to the HFS of the fip 
coincide: 



'-^-^iteryP+str 



-2<V''xG^x AV^P >^^^= -2 < V^p xG^x AVf/' >^^^= (42) 



P(l-t;2) + l' 



Numerical value of the proton structure and vacuum polarization effects of the 27 amplitudes 



AE- 



HFS 



+ 2AE, 



HFS 



-0.021 meV. 



(43) 



yp,str "T ^'-^■^iteryP+str 

Hadronic vacuum polarization contribution to the HFS of the ground state in the fip was 
studied in Ref. [3^. Here we present it in the different form using the expressions (38) and 
(41): 

^HFS _ _^F r Pi^)ds 



AEi 



dkVvpik). 



(44) 



Dividing the integration range over s on the intervals where the cross section of the e^e~ 
annihilation into hadrons (p(s) = a^{e~^e~ — > hadrons)/3sa^^) is known from the experi- 
ment I34II we can make numerical integration in Eq.(44). The result coincides with obtained 
in R.ef . l33l : 



AE 



HFS 
HVP 



0.004 meV. 



(45) 



11 




V. PROTON STRUCTURE EFFECTS, SELF ENERGY AND VERTEX 
CORRECTIONS OF ORDER a{Zaf 

There exists real number of important contributions of order which are presented in 
Fig. 4, 5. Radiative corrections of these amphtudes including recoil effects were studied earlier 
both in the Lamb shift and HFS of the hydrogen-like systems 3^, 3^. Radiative photons 
were taken in the Fried- Yennie (FY) gauge j37, 3^ 3^ where the mass shell amplitudes don't 
contain infrared divergences. Infrared finiteness of the Feynman diagrams in this gauge 
gives the possibility to make standard subtraction on the mass shell without introducing the 
photon mass. Let us consider radiative corrections which are determined by the self-energ 
insertions in the muon line. The renormalizable mass operator in the FY gauge is equal [U, 



= -{P - / dx -. (46) 

vr Jo mfx + (mf — p^)(l — x) 

Making the insertion (46) in the lepton tensor of the two-photon exchange diagrams and 
using the projection operators (35) we can construct the hyperfine splitting part of the 
quasipotential for the amplitudes in Fig. 4. In this case as before the vertex of the proton- 
photon interaction is determined by electric and magnetic form factors because the typical 
loop momenta are of order the muon mass. The contraction of the lepton and proton tensors 
over the Lorentz indices and the Dirac 7 matrix trace calculation were made in the system 
Form [4^. In the Euclidean space of the variable k we can present the correction to the 
HFS of the muonic hydrogen as follows: 

^K^iE = ^^"^^ '^^0- /' xdx r kdk r sin^ (t)dWsE{k, 0, x), (47) 

" TT^n'^ TT Jo Jo Jo 

' (fc^ -|- 4m2 cos^ (/))[(xmf + x/c^)^ + 4mfx^/c^ cos^ 0] ^ ^ 
^k'^F^ix + 6x) cos^ ^A;^xF| + 16m?F2 cos^ (piAFix - F2X - 2F2x) + 

m2 1712 

+ lQml cos^ (j){F^x + 6F^x + AF1F2X + 8F1F2X + F|x + 2F|x) + 
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abed 

FIG. 5: Proton structure and muon vertex effects of order a{Zctf' . 



+32m?a;Fi(Fi + F^) ^Fl^cos^ 



^Flx - l^k^Flx cos^ 0+ 



+16A;2xcos2 + 4F1F2 + F|) + 'i1k^F^x{F^ + F2)|. 

After analytical integration over the angle we present the contribution (47) in integral 
form which was used for numerical calculation: 



1 rao 

(5;o / xdx I dk< 







-^ + 32Fi(F, + F2) 
mo 



1 



+ 

2km\ 



k^Fi 



Qmjk^F,^ 



m 



4A; 
mi 



hi{k^ x) 
k ^ 



(49) 



2 F2(2Fi + F2)x f(xmf + xr) 

m2 m2 



\h2{k,x) hi{k,x) 
e ( 1 



+ 



k 



h2{k,x) mKxmf + xfc^) \h2{k,x) hi{k,x) 
hi{k, x) = kJ Amlx'^k'^ + (xmf + xk'^y + (xm^ + xfc^) \/4m2 + A;^, 



/i2(/i;, a;) = y^4mfx^A;^ + (xmf + xA;^)^ + (xm^ + xA;^). 
Numerical value is equal 

AE^^% = 0.008 meV. (50) 

Let us consider calculation of the vertex corrections. The renormalizable expression of 
the one-particle vertex operator in the FY gauge was obtained in Ref. ji^l (Pi = ml): 



A^ip,p-k) 



a 



»1 rl 

dx dz 
A-n Jo Jo 



M _|_ M 



A 



A2 



(51) 



where A = mfx + 2pk{l — x)z — k'^z{l — xz), the functions F^p ^ Fjf^ were determined in 



Ref. 



MB 



41|. The lepton tensor can be divided into two parts: 

v{pi)Fl;^\—pi — k + mi)7^t>(gi)(A;2 — 2/c°mi)[m^x — k^z{l — xz) + 2miA; 
{k'^ — 4A;o^f)[(mfx — k'^z{l — xz))^ — ArriikQx'^ z"^] 



^x^l 



(52) 
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^(0(2) ^ v{Pi)Fj)''\-pi -k + mi)^^v{qi){k^ - 2A;°mi)Ka; - ez{l - xz) + 2mik^x^]'' 

(k^ - Akimj) [(mix - k^z{l - xz)f - Amlklx'^z'^Y 

(53) 

Remaining for the simplicity the main contribution over mi/m2 we write this type vertex 
corrections as follows: 

^ Vi(j^.fcj.)|Fi(Fi + f2)-(l + K)l 

(/c^ + 4mi cos^ 0)(A:^ + Am^ cos^ 0) [[mfa; + kP'z{l — zx)Y + 'Arnlk'^ cos^ 0x^2;^] ' 

k, 0) = -2m^a;2(l - + k^ml{Qx^z^ - %x^z^ - ^x^z + ^xz - 3x) + (55) 
+k^{Ax^z^ - Qx^z^ - hx^z^ + 12x^3 _ 2a;^2 _ g^2 ^ 3^-)^ 

. - (Iff^ I .(1 - .)<te I . W /; . (56) 



.TT/ (1 + k)?™^ io 

F2(a;,fc,(/))Fi(Fi + F2) 



{k^ + 4mf cos^ 0)(A;^ + Aw^ cos^ 0) [[mfx + k'^z{\ — zx^Y + Amfk'^ cos^ 0a;^2;^]^ ' 

1^2(2;, k, 0) = mia;^^(2^ - 1) - k'^mjxz'^iAxz'^ - 2xz -Az + 2)+ (57) 
+k'^z^{2x'^z^ - x^z^ - Axz"^ + 2xz + 2^-1). 
The iteration contribution is equal 

A<ef 2, .ert =< V,, X X V,, >fjf = F^^ (^)7^' ^ dx f dk^ , (58) 

After analytical integration in Eqs.(54) and (56) over the angle and subtraction (58) (one 
photon is the Coulomb-like and the other one contains the hypcrfine part of the potential 
with the value of magnetic form factor at zero point) we have the expressions of the diagrams 
(a) and (b) in Fig. 5: 



Jo \8k{l + K)mim2x'^z'^ 

X {Gx^z"^ - Sxz^ - ?,xz + 8^ - 3) + k^z{4:X^z^ - 6x^z^ - bx^z^ + Uxz"^ - 2xz - 6z - 3) 

Fi{Fi + F2)x 



(59) 

X 



+ 



c 



2 



+ 



2(1 + K)mim2kx^z' 



4 



m\x'^z{2z - 1) - 2k'^mlxz^{2xz^ - xz - 2z + 1) + k'^z^{2x^z^ - x'^z^ - Axz"^ + 2xz + 2^-1) 



X 



X 



6(a2-62)(ft2_c2)2 a(a2-62)(a2_c2) 2cy/lT^{b' - c'){o? - c^) 

4/il 



v^TT? a/TTc^ ^/^T 
J h 

I /70 o\/o o\i 



2c3(62_c2)(a2-c2) c(62 -c2)(a2 -c2) c(62-c2)(a2-c' 



2^2 



^ F f 
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,2 



a 



fc^ 2 [m\x + k'^z{l — xz)] 



= ' ' .^2.L.., " ■ (60) 



4mf' 4m2' Amlk'^x'^z'^ 

Numerical value of vertex correction (59) is equal 

Ai^f/Li = -0.014 meV (61) 

Next vertex type diagram with one rounded photon and two exchanged photons is the 
diagram of the "jellyfish" type. Its contribution to the energy spectrum is of order a{Za)^. 
At small loop momenta this diagram gives the finite answer in the FY gau ge. The lepton 
tensor relating to the diagrams (c) and (d) in Fig. 5 was obtained in Ref. [36|: 



n=l 



where A has the form as in Eq.(51). The tensor functions are written explicitly 

in Ref. 3^|. The character of further transformations of the amplitudes (c), (d) in Fig. 5 



to construct the HFS part of the potential is the same as for other amplitudes shown in 
Fig. 4, 5. Omitting the details of such transformations which were carried out by means of 
analytical system Form [40] we write here three contributions to the HFS corresponding to 
the functions M^"^ in the leading order over mi/m^. 



A<^f.,/... = f\dx f\l-z){l-3xz) rkdkF,{F, + F,) x (63) 

' TT'^ri'^ JO Jo Jo 

sin^ (f)d(f) [mix + fc^z(l — xz)] 



{k"^ + 4m2 cos^ 0) [mf X + k'^z{l — xz)]'^ + Arriik'^ cos^ (px'^z'^ ' 

Ai?S,,.. = r ,dx [\l - z)dz r kdkF.iF, + F.) X (64) 

■' 6TX-^n'^ Jo Jo Jo 



sin^ [m\x + k'^z{l — xz)]'^[k^xz^{l — xz) + m\[x'^z + 2xz — x — ?>z)] 

{k^ + Ami cos^ (j)) { [m\x + k'^z{l - xz)]'^ + Amfk^ cos^ (px^z^^ 

'^ES%,„., = i'd - ^)z'd^ f k'^dk'n!F,iF, + F.) X (66) 

2 \ r sin^ (pdcj) [m\x + k'^z{l — xz)Y 

' Jo (A;2 + 4micosV){[m?x + A:2z(l-xz)]2 + 4mffc2cos2 0x222}3 

The integration over the angle can be done in Eqs.(63)-(65) analytically. Omitting inter- 
mediate expressions we can write final numerical result to the HFS of the /ip: 

A^S-/^ = E ^E^,7euyf^sh = 0.004 meV. (66) 

n=l 

In the point-like proton approximation when the nucleus form factors entering the Feynman 
amplitudes in Fig. 4,5 are changed on their values at /c^ = (-Fi(O) = 1, F2{0) = n) the 
contributions (63)-(65) will increase twofold. 
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TABLE I: Corrections of orders a^, to the ground state HFS in the muonic hydrogen. 



Contribution to HFS of fip 


Numerical value in meV 


Reference 


The Fermi energy 


182.443 


[18j, (12) 


Muon AMM correction a^E^ of order a^,a^ 


0.213 


fl8] 


Relativistic correction ^{Za)^E^ of order a® 


0.015 


[43] 


Relativistic and radiative recoil corrections 
with the account proton AMM of order 


0.014 


[30] 


One-loop electron vacuum polarization 
contribution of I7 interaction of order 


0.374 


(18) 


One-loop muon vacuum polarization 
contribution of I7 interaction of order 


002 


(19) 


Vacuum polarization corrections of orders a^, 
in the second order of perturbative series 


0.736 


(30)+(33) 


Proton structure corrections of order 


-1.215 


[22j, (40) 


Proton structure corrections of order 


-0.014 


[8] 


Electron vacuum polarization contribution-l- 
proton structure corrections of order 


-0.021 


(43) 


Two-loop electron vacuum polarization 
contribution of I7 interaction of order 


0.003 


(21)+(24) 


Muon self energy -|- proton structure 
correction of order 


0.008 


(50) 


Vertex corrections -|- proton structure 
corrections oi order a 


-U.Ui4 


(blj 


"Jellyfish" diagram correction -|- 
proton structure corrections of order 


0.004 


(66) 


HVP contribution of order 


0.004 


(45) 


Proton polarizability contribution of order 


0.084 


[16] 


Weak interaction contribution 


0.002 


[44] 


Summary contribution 


182.638 ± 0.062 





VI. CONCLUSION 

We made the calculation of different quantum electrodynamical effects, effects of the pro- 
ton structure and polarizability, the hadron vacuum polarization to HFS of muonic hydrogen. 
The corrections of order and were considered. Working with the vacuum polarization 
diagrams we take into account that the ratio /la / mg is very close to 1 and don't increase the 
order of corresponding contributions. Obtained numerical results are presented in the Table 
1. We include here also QED corrections to the Fermi energy which are determined by muon 
anomalous magnetic moment a^E"^ uA (experimental value of muon anomalous magnetic 
moment a^^^ = 11659203(8) x 10~^° |42l was used), the Breit relativistic correction of order 
{Za)^ ji^], relativistic and radiative recoil effects of the same order [Za)^mi/m2 with the 



16 



account of the proton anomalous magnetic moment j30|, the proton structure corrections of 
order {Zaf\a.{Zaf' 0, the hadron vacuum polarization contribution js^ and the proton 
polarizability correction the weak interaction contribution due to Z boson exchange 

Let us point out some peculiarities of this investigation. 

1. The effects of the vacuum polarization play very important role in the case of the 
muonic hydrogen. They lead to essential modification of the spin-dependent part of the 
quasipotential of the one-photon interaction. 

2. We took into account consistently the proton structure in the loop amplitudes by 
means of electromagnetic form factors. The point-like proton approximation gives essentially 
increased results (approximately twofold). 

3. The calculation of muon self-energy and vertex corrections of order aiZaf was done on 
the basis of the expressions for the lepton factors in the amplitude terms of the quasipotential 
obtained by Eides, Grotch and Shelyuto. We supplemented these relations by the subtraction 
of the iteration terms of the potential. 

Total value of the ground state HPS in the muonic hydrogen shown in the Table 1 can 
be considered as definite guide for the future experiment which is prepared [23]. Numerical 
values of the corrections were obtained with the accuracy 0.001 meV. Theoretical error 
connected with the uncertainties of fundamental physical constants (fine structure constant, 
the proton magnetic moment etc.) entering the Fermi energy compose the value near 10~^ 
meV. Other source of theoretical uncertainty is connected with the corrections of higher 
order. Its estimation can be found from the leading correction of the next order on a and 
mi/m2 in the form: a(Za)^ln(Za)^/7r ^ 0.0005 meV ( the value of fine structure constant 
is = 137.03599976(50) Q). 

It is useful to compare the summary result for the ground state HFS in the muonic 
hydrogen obtained in this work (see the Table 1) with that one which can be founded in the 
point like proton approximation when we take into account only the values of the proton 
electromagnetic form factors at fc^ = 0: Ge(0) = 1, G'a/(0) = /i (with the exception of 
the Zemach correction). In this approximation the result of the ground state HFS may be 
presented with the accuracy 0((mi/m2)a^) as follows [l8l |: 

(67) 



AE'^f'iQED) = E^jl - 2/xai?p + ^{Zaf + a^ + a{Za) (\n2 - ^) + 



1 



1 + K 



— _2 I2 In 3^ + {Zaf^^ ( 2(1 + k) + — \n{ZaY^- 



TT ml — ml mi ' ' mim2 



.(12-11.) ^311 ^.(11 + 31.) 



4 / 18 36 

181.177 meV. 



-'^{Zaf\n{Za)-^mlrl 



Essential difference between this numerical value and 182.638 meV obtained in our study 
can be explained by some reasons: the modification of the Breit potential due to electron 
vacuum polarization for the muonic hydrogen, the effects of the proton structure in the 
two-photon and three-photon interactions, hadronic vacuum polarization and proton polar- 
izability effects in our calculations. Further improvement of theoretical result presented in 
the Table 1 is connected first of all with the corrections on the proton structure and po- 
larizability which give the theoretical error near 340 ppm. The most part of this error is 
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determined by the proton structure corrections of order iZaf' (the Zemach correction). So, 
the measurement of the hyperfine sphtting of the levels IS* and 25* in the muonic hydrogen 
with the accuracy 30 ppm will lead to more accurate value (with relative error 10~^) for 
the Zemach radius which than can be used for the improvement of theoretical result for 
the ground state hydrogen hyperfine structure and more reliable estimation of the proton 
polarizability effect. The increase the number of the tasks due to excited states of simple 
atoms and the inclusion new simple atoms where the hyperfine structure of the energy 
spectrum is studied will decrease the uncertainties in the determination of physical funda- 
mental parameters and increase the accuracy for the check of the Standard Model in low 
energy physics. 
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